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We have studied the electronic properties in aperiodic graphene superlattices of Thue-Morse 
sequence. Although the structure is aperiodic, an unusual Dirac point (DP) does exist and its 
location is exactly at the position of the zero-averaged wave number (zero-fc). Furthermore, the 
zero-fc gap associated with the DP is robust against the lattice constants and the incident angles, 
and multi-DPs can appear under the suitable conditions. A resultant controllability of electron 
transport in Thue-Morse sequence is predicted, which may facilitate the development of many 
graphene-based electronics. 



Graphene has attracted enormous attention of experi- 
mentalists and theorists [ll-Q since its discovery. The in- 
terest is driven by its potential technological applications 
and unconventional low energy behavior. This is mainly 
because graphene has a unique band structure with 
the conductance and valance bands touching at Dirac 
point (DP). Recently, scientists anticipate that graphenc- 
based electronics may supplement silicon-based technol- 
ogy, which is nearing its limits For superlattices are 
vastly successful to control the electronic transport [lol|. 
to facilitate the application of graphene-based electron- 
ics, the graphene superlattices (GSLs) with electrostatic 
potential or magnetic barrier have also received broad fo- 
cuses of theoretical and experimental investigations [lll- 
[iol . In such GSLs, the new DP appears in the band struc- 
3] and it's exactly located at the energy cor- 
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responding to zero-averaged wave number (zero-fc) 
Contrary to Bragg gaps, the zero-fc gap associated with 
the new DP is insensitive to both the lattice constant and 
the structural disorder, resulting in better controllability 
of electronic transport in GSLs. Most recently, Zhao and 
Chen predicted a controllable electron transport in a Fi- 
boncacci quasi-periodic GSl[i§]. 

In this letter, we investigate the electronic band gaps 
and transport in the GSLs of Thue-Morse (TM) se- 
quence. Although such a TM sequence is aperiodic and it 
is more "disordered" than the Fibonacci sequence [illli^ . 
we find that zero-fc gap and DP do exist in such a GSL 
TM sequence, which results in robust electronic trans- 
port properties. 

For a TM sequence, it can be generated in various 
ways. For an n-th TM sequence, S{n), it contains 
only 2 elements A and B, and follows the inflation rule: 
A — >■ AB, B — > BA from generation to generation with 
5(1) = AB (for example, see Ref. [H]). Let S^^^{n) and 

'^i/2('"') denote the left and right half parts of an n-th 
order TM sequence respectively, then iteration relation 



for any {n + l)-th TM sequence is written as 

S{n + l) = S{n) S^/,{n) S^/^{n). (1) 

Naturally, 5(2) = ABBA, 5(3) = ABBABAAB, and so 
on. In our cases, A [B] denotes barrier Va {Vb) with its 
width wa (wb)- It is readily to know that the number 
of A, Na, equals to that of B, Nb, in any S{n), i.e., 
Na = Nb = 2"-\ 

Meanwhile the charge carriers near the K point in 
graphene (near the Fermi level) are described by its 
Hamilton: H = vpa ■ p -|- V{x)I, where V{x) is a barrier 
or well, vp PS lO^m/s is the Fermi velocity, a = {ax,^y) 
are Pauli matrices, and / is a 2 x 2 unit matrix. The 
solution of H acting on the electronic pseudospin wave- 
functions leads to a transfer matrix 17 



/ cos{qjAx-ej) ^sin(qjA£) \ 
Mj{Ax,E,ky)= I . sh°{ljAx) cos{q°Ax+0.i} I 1 (2) 



which connects the wave functions at two positions x and 
X -\- Aa; inside the jth potential with 9j = aicsiniky/kj), 
here ky is the y component of the wave number kj, 

and qj = sign(fcj), /fc^ — /c^ for k'l < k^ , otherwise 



= k} — k'i for fc.^ > fc?. The electronic transmis- 



qj 1. Y '-y '-J ^"^ ■•'y - '"J ■ 

sion coefficient t{E, ky) in such devices can be obtained 

by [i3 



t{E,ky) 



2 cos ^0 



(x22e -I- xiie*^^) - a;i2e' 



X2l 



where 9o{9e) is the incident (exit) angle, and Xij(i,j 



(3) 



N 



1, 2) is the element of X[5(n)] = JJ Mj{wj,E,ky), which 

is the entire transfer matrix of a GSL TM sequence. If we 
let Xn = Tr{X[5(7i)]}, it is readily to derive the iteration 
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FIG. 1: (Color online) Trace-maps for GSL TM sequences of 
variable order n under (a) wa,b = 20 nm, (b) wa,b ~ 25 nm 
and (c) wa ~ 2wb = 40 nm with Va = 50 meV, Vb ~ 
and the wave number ky = 0.015 nni^^. (d), (e) and (f) are 
band structures for n — 8 corresponding to (a), (b) and (c) 
respectively as ky varies. The horizontal dashed black line 
denotes the new Dirac point's location, and the white areas 
are band gaps and the solid areas are passing bands. 



relation for the trace map of the n-th GSL TM sequence 
as follows [ii-il: 



^1-2) 



(4) 



Treating an n-th TM sequence as a unit cell, from 
Bloch's theorem, we have cos(/3a;A„) = 2;„/2, where 
An — Nawa + Nbwb- From Eq.Q, we can calculate 
the change of the trace map as a function of the order n. 

In Fig. [T] (a) and (b), we plot the trace maps for two 
kinds of GSL TM sequences dependent on the order n, 
at the incident angle = 10°. We take wa,b = 20nm in 
Fig. 1(a), and wa,b = 25nm in Fig. 1(b). We find that 
the passing bands are split into more and more sub-bands 
as n increases; when n > 6, the band structures almost 
become the discontinuous bands. However, we note that 
the center position of the Gap A in Figs. 1(a) and 1(b) 
are the same, and the positions of other gaps are shifted 
with the change of the lattice constant. Corresponding 
to Fig. 1(a) and 1(b), the electronic band structures for 
a TM sequence of n = 8 have been shown in Fig. 1(d) 
and 1(e), respectively. It is seen that a new DP appears 
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FIG. 2; (Color online) Transmission spectrum for an 8-th 
GSL TM sequence (256 layers) where Va = 50 meV, Vb = 
0, {a.)wA ~ 15 nm (solid blue line), 20 nm (dashed green), 
25 nm (dotted red) with a fixed ratio wa/wb ~ 1 and an 
incident angle Oq — 10°, and (b) ~ 10° (solid blue line), 
20° (dashed green), and 25° (dotted red) with fixed lattice 
constants wa,b = 20 nm. 



inside the Gap and this new DP actually locates at 



N 



N 



(5) 



This condition is valid for periodic and aperiodic 
graphene superlattices Since Na — Nb in our 

cases, it is easily to obtain the energy for fc = as follows: 



E 



Va + Vb ■ wb/wa 
I + wb/wa 



(6) 



In the cases of Fig. l(a,b,d,e), it is given hy E = 25meV 
as WA = wb- 

Since the center of the Gap A in Fig. 1 is located at 
zero-fc, which is denoted by a dashed black line in FigH] 
(a) and (b), we may call it as the zero- A: gap. According 
to Eq.(jni), the position of the zero-fc gap for the GSL TM 
sequences is not shifted with the lattice constant itself, 
but is shifted with the change of the ratio of wa/wb- As 
shown in Fig [I](c) and (f), the position of the zero- A: gap 
and the DP move to « 33.33meV for wa/wb = 2. 

FiglH (a) shows the effect of the lattice constants' vari- 
ation with a fixed ratio of wa/wb = 1- Based on Fig [2] 
(a), it is obvious that the zero-fc gap is insensitive to the 
lattice constants themselves. However, the positions of 
other gaps with higher energy are highly dependent on 
the lattice parameters. From Fig. [H^b), one can also 
find that the position of zero-fc* gap is weakly dependent 
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FIG. 3: (Color online) Electronic band structures for (a) 
WA,B ~ 20 nm, (b) ■wa,b = 40 nm and (c) wa,b = 55 nm; 
(d) band gaps depending on the lattice constant w with fixed 
0.015 nm~^ and wa = wb = w. Other parameters are 



ky 

the same as those in Fig.l. 



on the incident angle 6*0 while other Bragg gaps change 
sensitively with Oq. 

Furthermore, one has known that the multi-Dirac- 
points could appear in the GSLs with periodic poten- 
tial structures [16, 
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Here we point out that the extra 
Dirac points, located at ky ^ 0, could also emerge in the 
GSL TM sequence as the lattice constant increases. See 
FigslSlJa) and (b), the slope of the band edges near the 
DP gradually turns smaller as wa (wb) increases from 
20 nm to 40 nm. When wa (wb) is larger than 40 nm, 
as that shown in Fig. [3] (c), the additional Dirac points 
appear. From FigsISl^a) to (c), we may also find the ap- 
pearance of additional Dirac points associating with the 
variation of the zero-fc gap. In Fig (3] (d), we show the 
energy band gaps as a function of the lattice constants 
in the case of wa = wb- In the process of the appear- 
ance of the additional Dirac points, the zero-fc gap opens 
and closes oscillationally while the other gaps are shifted 
greatly. 

Finally, we calculate the total conductance G 
[H and the Fano factor F ^ in GSL TM se- 



quences, which are given by G=Gq /J 

tt/2 



IT 12 



TcosOqcIBq and 



F 



=j:i%T{l^T)cos9od0o/j:i^-/, 



T cos9od6o, where 



T=|f| and Go='2e mvpLy/h , with Ly denoting the 
width of the graphene stripe in the y direction. We 
present the total conductance and the Fano factor 
as a function of Fermi energy with different lattice 
constants in FigHl What we should notice is that the 
angular-averaged conductance G's curve reaches its 
minimum at the DP and forms a linear cone around the 
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FIG. 4: (Color online) Conductance ( (a) and (c) ) and Fano 
factor ( (b) and (d) ) vs Fermi energy in an 8-th GSL TM 
sequence. Other parameters are the same as those in Fig.l. 



DP, and F at the DP's location reaches the value of 
1/3 approximately [20I. \2l\. Moreover, the conductance 
and the Fano factor shift with the ratio of wa/wb since 
the DP's location relies on this ratio. Therefore, this 
suggests that the conductance of GSL TM sequence 
could be modulated by the ratio of the lattice constants. 

In summary, we have studied the electronic trans- 
port properties in the GSLs of Thue-Morse aperiodic se- 
quence. It is found that an unusual DP does exist and 
its location is exactly at the energy corresponding to the 
zero-fc in such aperiodic structures. The zero-fc gap asso- 
ciated with the DP is robust against the lattice constants 
and the incident angles, and the extra-DPs can appear 
under the suitable conditions. A resultant controllable 
electron transport is predicted, and our results should 
facilitate the development of graphene-based electronics. 
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